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ABSTRACT 

Planck is expected to answer crucial questions on the early Universe, but it also provides further understanding 
on anomalous microwave emission. Electric dipole emission from spinning dust grains continues to be the favored 
interpretation of anomalous microwave emission. In this paper, we present a method to calculate the rotational 
emission from small grains of irregular shape with moments of inertia Ii > I2 ^ la- We show that a torque-free 
rotating irregular grain with a given angular momentum radiates at multiple frequency modes. The resulting spinning 
dust spectrum has peak frequency and emissivity increasing with the degree of grain shape irregularity, which is 
defined by Ii : I2 '■ la- We discuss how the orientation of dipole moment fi in body coordinates affects the spinning 
dust spectrum for different regimes of internal thermal fluctuations. We show that the spinning dust emissivity for the 
case of strong thermal fluctuations is less sensitive to the orientation of fi than in the case of weak thermal fluctuations. 
We calculate spinning dust spectra for a range of gas density and dipole moment. The effect of compressible turbulence 
on spinning dust emission is investigated. We show that the emission in a turbulent medium increases by a factor 
from 1.2-1.4 relative to that in a uniform medium, as sonic Mach number Ms increases from 2-7. Finally, spinning 
dust parameters are constrained by fitting our improved model to five-year Wilkinson Microwave Anisotropy Probe 
cross-correlation foreground spectra, for both the Ha-correlated and 100 /xm-correlated emission spectra. 

Subject headings: ISM: dust, extinction — ISM: microwave emission — galaxies: ISM — infrared: 
galaxies 



1. INTRODUCTION 

Cosmic Microwave Background (CMB) experiments 
(see Bouchet et al. 1999; Tegmark et al. 2000; Efs- 
tathiou 2003; Bennet et al. 2003) are of great importance 
for studying the early universe and its subsequent ex- 
pansion. Precision cosmology with Wilkinson Microwave 
Anisotropy Probe ( WMAP) and Planck satellite requires 
a good model of the microwave foreground emission to 
allow the subtraction of microwave foreground contami- 
nation from the CMB radiation. 

Three well-known components of the diffuse microwave 
Galactic foreground consist of synchrotron emission, free- 
free emission from plasma (thermal bremsstrahlung) and 
thermal emission from dust. However, in the microwave 
range of frequency from 10-100 GHz, anomalous emis- 
sion which was difficult to reconcile with the components 
above was reported (Kogut et al. 1996a, 1996b). 

An explanation for the anomalous emission was pro- 
posed by Draine & Lazarian (1998ab, hereafter DL98 
model), where it was identified as electric dipole emis- 
sion from rapidly spinning dust grains (hereafter spin- 
ning dust emission). Although spinning dust emission 
had been discussed previously (see Erickson 1957; Fer- 
rara & Dettmar 1994), Draine & Lazarian were the first 
to include the variety of excitation and damping pro- 
cesses that are relevant for very small grains. 

While the DL98 model appears to be in general agree- 
ment with observations (see Lazarian & Finkbeiner 2003; 
Finkbeiner 2004), it did not account for a number 
of physical effects, namely, the non-sphericity of grain 
shapes, internal relaxation, and transient spin- up due to 

^ Astronomy Department, University of Wisconsin, Madison, WI 
53706, USA 

^ Department of Astrophysical Sciences, Princeton University, 
Princeton, NJ 08544, USA 



ion collisions. 

The Planck Collaboration 2011b has reported new ob- 
servations of spinning dust emission in new environ- 
ments. Spinning dust now provides a potential diag- 
nostic tool for interstellar dust properties (shape, size 
distribution, dipole moment). A comprehensive model 
of spinning dust for different grain shapes is required. 

Recent studies showed that the correspondence of the 
DL98 model to observations can be improved by adjust- 
ing the parameters of the model. For instance, the five- 
year ( WMAP) data showed a broad bump at ^ 40 GHz in 
the Ha-correlated emission (Dobler & Finkbeiner 2008; 
Dobler, Draine & Finkbeiner 2009). The - 40 GHz fre- 
quency of the peak is higher than the value predicted 
(^ 23 GHz) by the DL98 model for standard parame- 
ters of the warm ionized medium (WIM). Dobler et al. 
(2009) showed that the bump is consistent with the DL98 
model modified so that the characteristic dipole moment 
of grains is decreased and gas density of the WIM is in- 
creased, relative to the typical spinning dust parameters 
in Draine & Lazarian (1998b, hereafter DL98b). 

Ali-Hai'moud, Hirata & Dickinson (2009) revisited the 
spinning dust model and presented an analytic solution 
to the Fokker-Planck (FP) equation that describes the 
rotational excitation of a spherical grain if the discrete 
nature of the impulses can be neglected. 

Hoang, Draine & Lazarian (2010, hereafter HDLIO) 
improved the DL98 model by accounting for a number 
of physical effects. The main modifications in their im- 
proved model of spinning dust emission are as follows. 

(i) Disk-like grains rotate with the grain symmetry axis 
ai not perfectly aligned with the angular momentum J. 
The disaligned rotation of a; with ai causes wobbling of 
the grain principal axes with respect to w and J due to 
internal thermal fluctuations. 
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(ii) Distribution functions for grain angular momen- 
tum and velocity are obtained exactly using the Langevin 
equation (LE) for the evolution of grain angular momen- 
tum in an inertial coordinate system. 

(iii) The limiting cases of fast internal relaxation and 
no internal relaxation are both considered for calculation 
of the angular momentum distributions and emissivity of 
spinning dust. 

(iv) Infrequent collisions of single ions which deposit 
an angular momentum larger than the grain angular 
momentum prior to collision are treated as Poisson- 
distributed events. 

The wobbling disk-like grain has anisotropic rotational 
damping and excitation. Such an anisotropy can increase 
the peak emissivity by a factor ~ 2, and increases the 
peak frequency by a factor 1.4 — 1.8, compared to the 
results from the DL98 model. 

The effect of the grain wobbling on electric dipole emis- 
sion was independently studied in Silsbee et al. (2011) 
using a Fokker-Planck (FP) equation approach, but they 
disregarded the effect of internal relaxation and transient 
spin- up by infrequent ion collisions. 

Earlier models of spinning dust emission dealt with ax- 
isymmetric grains having moments of inertia Ii > I2 = 
I3 (see HDLIO; Silsbee et al. 2011). An axisymmetric 
grain of an angular momentum J radiates in general at 
four frequency modes 4) ± t/j.ii) and ■;/', where (f> = J/I2 
is the angular frequency of precession of the grain sym- 
metry axis ai about J, and ip = J cos 6(1/ Ii — I//2) 
with being the angle between J and ai is the angular 
frequency of the grain rotation about its symmetry axis 
ai. For an "irregular" (i.e. triaxial body) grain with 
Ii > I2 > h, in addition to the precession of the grain 
axis of maximum moment of inertia ai (hereafter axis of 
major inertia) or of the grain axis of minimum moment 
of inertia as (hereafter axis of minor inertia) around the 
angular momentum J, the grain principal axes wobble 
with respect to J. This wobbling, that is different from 
the wobbling due to the internal thermal fliictuations 
studied in HDLIO for the disk-like grain, can result in a 
more complex electric dipole emission spectrum, so that 
it modifies the spinning dust emission. We will quantify 
this effect in §3. 

Very small grains (smaller than 20 A) are important 
for spinning dust emission. But also within this range 
of grain size, the grain dust temperature has strong fluc- 
tuations due to absorption of UV photons (Greenberg 
1968; Draine & Anderson 1985). As a result, they can 
not be characterized by a single equilibrium temperature, 
and are described by a temperature distribution function 
(Draine & Li 2001). The temperature fluctuations induce 
wobbling of the grain axes with respect to the angular 
momentum and modify the spinning dust spectrum. This 
issue is addressed in §5.1. 

We investigate also spinning dust emission in the pres- 
ence of compressible turbulence. An increase of the spin- 
ning dust emission arises as a result of the non-linear 
dependence of emissivity on gas density. 

The seven-year WMAP data (Gold et al. 2011) and 
Planck early results (Planck Collaboration 2011ab) ap- 
pear to support spinning dust emission as a source of 
the anomalous microwave emission in new regions with 
a wide range of physical parameters. The WMAP ther- 



mal dust-correlated spectra in Gold et al. (2011) for two 
regions within the Galactic plane indicate that spinning 
dust emission should peak around 22 GHz. However, the 
improved model of spinning dust developed by HDLIO 
predicted a peak frequency above ~ 30 GHz, for all me- 
dia using the same physical parameters as in DL98b. We 
explore parameter space for the HDLIO spinning dust 
model, and characterize spinning dust parameters by fit- 
ting to the latest observation data. 

The structure of the paper is as follows. In §2 we 
present the assumptions and notations adopted through- 
out this paper. In §3, we describe the torque-free motion 
of an irregiilar grain and the exchange of vibrational- 
rotational energy. The; power spectrum and spinning 
dust emissivity for irregular grains are presented in §4. 
In §5, we investigate the influence of grain temperature 
fluctuations, dipole moment orientation, and explore the 
parameter space of gas density and magnitude of dipole 
moment for spinning dust emissivity. §6 is devoted to the 
effects of fluctuations of gas density due to compressible 
turbulence on spinning dust emission. Constraints on 
spinning dust parameters obtained by fltting theoretical 
model to the Ha-correlated and thermal-dust-correlated 
emission spectra are presented in §7. Discussion and 
summary of principal results are given in §8 and 9, re- 
spectively. 

2. ASSUMPTIONS AND NOTATIONS 

2.1. Assumptions 

In the present paper, we follow the same assumptions 
for grain size distribution and rotation dynamics as in 
HDLIO. The grain size a is defined as the radius of a 
sphere of equivalent volume.'^ Unless stated otherwise, 
we assume throughoiit the paper that grains smaller than 
02 = 6 X 10~^ cm are planar, and grains larger than 02 
are approximately spherical. We adopt rotational damp- 
ing and excitation coefficients due to gas-grain interac- 
tion, infrared emission, and electric dipole damping from 
HDLIO. Our notation is summarized in Table 1. 

2.2. Grain electric dipole moment 

The electric dipole moment of a grain arises from the 
intrinsic dipole moment of molecules within the grain and 
from the asymmetric distribution of the grain charge. 
The former is shown to be dominant (see DL98b). 

The grain dipole moment /x can be written as 



(1) 



where /x^ is the intrinsic dipole moment of an uncharged 
grain, Ze is the grain charge, ax is the excitation- 
equivalent radius of the grain, and the vector ea^ is the 
displacement between the grain center of mass and the 
charge centroid. 

Following DL98b, the magnitude of the dipole moment 
is given by 



/i^ = 23 



(^)'(Z2) + 3.8 



0.4 D 



a-7 



-7 D^ (2) 



where {Z"^) is the mean square grain charge, j3 is the 

* Denote V be the volume of the grain, then a is defined as 
V = 4na^/3. 
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dipole moment per atom of the grain, and a_7 
a/10~^cm. Above, e ~ 0.1 is assumed. 



TABLE 1 
Notations and Meanings 



Symbol 



Meaning 



gram size 
frequency 
H nucleus density 
H ionization fraction 
metal ionization fraction 
gas temperature 
dust vibrational temperature 
decoupling temperature for 
vibrational-rotational energy exchange 
gas damping time 
electric dipole damping time 
characteristic dipole moment 
electric dipole moment 
thermal angular velocity 
angular momentum 
basis vectors of grain principal axes 
eigenvalues of moment of inertia tensor, 
coordinate systems associated to J||z 
inertial coordinate system fixed to the lab 
Euler angles of grain axes in external 
system xyz 
2IiErot/ with -Erot rotational energy 
emission power 
emissivity at frequency u 
oriented so that /ij = /ij = 
fl perpendicular to ai 
25% of grains having /3 = /9o/2, 
50% of grains having /3o 
and 25% grains having 2/3o 



V 

XH 

xm 

Tgas 
Tvib 
Td 

TH 

Ted 

/3o 

J 

ai, a2, as 
Ilj2,h 
X, y,z 

61,62,63 

e, (j>, i) 
1 

Pod 

case 1 
case 2 
model A 




Fig. 1. — Torque-free motion of a grain having three principal 
axes ai,a2 and aa in an inertial coordinate system with z||J and 
X, y ± J described by three Euler angles 9, (f> and xp. 

Because the orientation of in the grain body is not 
well constrained, in the present paper, we study two lim- 
iting cases. In case 1, is fixed in the grain body such 
that ^1 = ^2 = — In case 2, /x is perpendicular 

to the grain axis of major inertia ai, so that ^\ — and 
;2 = ^2^/2 ^ ^"^12, where = + Mi- 
Table 2 presents typical physical parameters for vari- 
ous phases of the interstellar medium (ISM), including 
the cold neutral medium (CNM), warm neutral medium 
(WNM) and warm ionized medium (WIM). 



TABLE 2 

Idealized Environments For Interstellar Matter 



Parameters 


CNM 


WNM 


WIM 


nj] (cm~^) 


30 


0.4 


0.1 


Tga= (K) 


100 


6000 


8000 


X 


1 


1 


1 


xn 


0.0012 


0.1 


0.99 




0.0003 


0.0003 


0.001 


y = 2n(H2)/nH 


0. 


0. 


0. 



The electric dipole moment, assumed to be fixed in 
the grain body system, can be decomposed into three 
components along the grain principal axes: 

/X = niai + fj.2a2 + Aisaa, (3) 

where ai , a.2 and a.3 are the basis vectors of principal axes 
of the grain. 



3. ROTATIONAL DYNAMICS OF GRAIN OF IRREGULAR 

SHAPE 

Very small dust grains (also PAHs) are expected to 
be planar, while larger grains are likely to be more-or- 
less spherical. Even though planar, small PAHs are not 
expected to be perfectly symmetric. In general, for both 
PAHs and larger particles, we expect the eigenvalues li of 
the moment of inertia tensor to be nondegenerate: Ii > 
I2 > ^3- We will refer to grains with nondegenerate 
eigenvalues as "irregular" . 

3.1. Torque-free rotation and flip states 

Let us consider an irregular grain characterized by a 
triaxial body having moments of inertia Ii > I2 > I3 
along the principal axes ai , a.2 and a^ . 

The dynamics of the irregular grain is more compli- 
cated than that of a disk-like grain. In addition to the 
precession of the axis of major inertia ai about J as in 
the disk-like grain, the axis ai wobbles rapidly, resulting 
in variations of the angle 9 between ai and J at a rate 9. 
Therefore, we expect the wobbling of the irregular grain 
can produce more frequency modes than those observed 
in the disk-like grain with I2 — I3 (sec HDLIO). 

To describe the torque-free motion of an irregular grain 
having a rotational energy i^rot, we use conserved quan- 
tities, including the angular momentum J, and a dimen- 
sionless parameter that characterizes the deviation of the 
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grain rotational energy from its minimum value, 



(4) 



Following Weingartner & Draine (2003), we define the 
total number of states s in phase space for q ranging from 
1 to g as 



s = 1 - 



where 



dip 



ipi = cos ^ 



h{h - h) + h{h - h) cos^ ip . 



1/2 



1/2 



(6) 



for q > Qsp and Vi = 7r/2 for q < qsp, with g^p = /1//2 
being the separatrix between the two regimes. 

The instantaneous orientation of the grain in an iner- 
tial coordinate system is given by three Euler angles ijj, <P 
and 9 (see Figure 1). The angular velocities are obtained 
by solving the Euler equations of motion. The torque- 
free motion of an irregular grain is treated in detail in 
Appendix A. 

Briefly, the Euler equations give rise to two sets of 
solutions corresponding to + and — signs for with 
i = 1,2 and 3 (see Eqs A4-A6). It can be seen that 
for q < q^p, two rotation states ± correspond to cji > 
and uji < 0, i.e., ai • J > and ai • J < 0. Following 
Weingartner & Draine (2003), we refer to these rotation 
states as the positive flip state and negative flip state 
with respect to ai (see also Hoang & Lazarian 2009a). 
For q > ^sp, similarly, two rotation states correspond to 
C1J3 > and W3 < 0, i.e., as • J < and as • J > 0. These 
rotation sates arc referred to as the positive and negative 
flip states with respect to as. 

3.2. Internal relaxation, thermal fluctuations and 
thermal flipping 

3.2.1. Internal relaxation and thermal fluctuations 

Internal relaxation (e.g. imperfect elasticity, Barnett 
relaxation; Purcell 1979; Lazarian & Efroimsky 1999) 
arises from the transfer of grain rotational energy to vi- 
brational modes. For cold grains, this process tends to 
result in nearly perfect alignment of the grain axis of 
major inertia with the angular momentum.^ Naturally, 
if the grain has nonzero vibrational energy, energy can 
also be transferred from the vibrational modes into grain 
rotational energy (Jones & Spitzcr 1969). 

For an isolated grain, a small amount of energy gained 
from the vibrational modes results in fluctuations of the 
rotational energy E'rot when the grain angular momen- 
tum is conserved. Such fluctuations in E^ot result in 
fluctuations of q and of the angle 6 between ai and J 
for an axisymmetric grain. For an irregular grain, the 
fluctuations in E^ot are described by fluctuations in s{q) 
(see Eq. 5). Over time, the fluctuations in Erot establish 
a local thermal equilibrium (LTE) at a rotational energy 
equilibrium temperature T^ot- 

* This rotation configuration has minimum rotational energy or 
highest entropy. 



3.2.2. Exchange of Vihrational-Rotational Energy 

The Intramolecular Vibrational-Rotational Energy 

Transfer process (IVRET) due to imperfect elasticity oc- 
curs on a timescalc 10^^ s, for a grain of a few Angstroms 
(Purcell 1979), which is shorter than the IR emission 
time. So, when the vibrational energy decreases due to 
IR emission, as long as the Vibrational-Rotational (V-R) 
energy exchange exists, interactions between vibrational 
and rotational systems maintain a thermal equilibrium, 
(5) i.e.. Trot ~ Tvib. As a result, the LTE distribution func- 
tion of rotational energy reads (hereafter VRE regime; 
see Lazarian & Roberge 1997): 



/vre(s,./) oc exp 



Etc 



B-trot 



exp 



Er, 



B-ivib 



•(7) 



The existence of V-R energy exchange depends on the 
rotational modes. In principle, the rotational energy 
can change in increments of huj, with lu being the an- 
gular frequency of rotational modes. However, at low 
vibrational energy Evih, the vibrational mode spectrum 
is sparse, and there may not be available transitions with 
Ai?vib ~ hurot- For the case ai parallel to J, the V-R en- 
ergy exchange occurs only when the grain is hot just after 
absorption of a UV photon, because the grain rotating 
with aill J has only one rotational mode w = Wrot = J/h- 
As the grain cools down further, there is no available vi- 
brational transition corresponding to an energy change 
hujy-ot, and the V-R coupling ceases. 

For disk-like or irregular grain, multiple rotational fre- 
quency modes are observed, and some modes correspond 
to energy huj > /kJrot (see Eq. 14 or Fig. 2). This 
may allow V-R energy exchange even when the vibra- 
tional energy is relatively low. Nevertheless, the energy 
separation between the vibrational ground state and the 
first vibrationally-excited state is large enough that V-R 
energy exchange is unlikely to take place after the PAH 
cools to the vibrational ground sate. The key, then, is 
to know the vibrational "temperature" of a cooling PAH 
at the time when V-R energy exchange ceases to be ef- 
fective. Sironi & Draine (2009) estimated a decoupling 
temperature Td = 65 K for a PAH with 200 C atoms. 

In regions with the average starlight background, the 
mean time interval between absorptions of starlight pho- 
ton will be ~ 10^ s for a PAH containing 200 C atoms. 
Because the PAH will cool below 65 K in ^ 10'^ s (see, 
e.g.. Fig. 9 of Draine & Li 2001), V-R energy exchange 
will be suppressed most of the time, taking place only 
during ^ 10^ s intervals following absorption of starlight 
photons. However, the 10^ s interval between photon ab- 
sorptions is short compared to the timescale over which 
J will change appreciably, and we will therefore approx- 
imate the V-R energy exchange process as continuous, 
even though it is episodic. 

Substituting ii^rot as a function of J and q from Equa- 
tion (4) into Equation (7), the distribution function for 
the rotational energy becomes 



/vre(s, J) 



Aexp y- 



qis)J' 
27ifcBTvib 



(8) 



where ^ is a normalization constant such that 
Jo /vre(s, J)ds = 1. 
Define a rotational temperature Tj = J^/2/ifcB. 
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Thus, the distribution for q becomes /vre « 
exp(-g(s)r,//Tvib). When Tj > Tvib, the grain is 
suprathermally rotating, and ai is ahgned with J. In 
contrast, when Tj <C Tvib, the grain rotates subther- 
mally, and thermal fluctuations randomize the grain axes 
with respect to J. For fixed angular momentum, the 
thermal fluctuations increase with Tyib- 

3.2.3. Thermal flipping 

Transfer of energy between the vibrational and rota- 
tional modes due to dissipation or thermal fluctuations 
will cause .Brot; therefore q, to vary with time. If q 
crosses the separatrix state q = qsp, the new flip state is 
effectively chosen randomly. Changes in q can also result 
from external forcing (e.g., collisions with gas atoms or 
absorption/emission of radiation), that may change both 
J and E^ot ■ If the external torques are not "chiral" , then 
the changes in q over time will result in equal probabil- 
ities of the two flip states. However, the chirality of the 
grain may lead to g'-changing processes that depend not 
only on q, but also on which flip state the grain is in. 
The present treatment neglects any grain chirality. 

4. ELECTRIC DIPOLE EMISSION FROM GRAINS OF 
IRREGULAR SHAPE 

Below we investigate electric dipole emission aris- 
ing from irregular grains. We first calculate emission 
power spectra arising from a torque-free rotating irregu- 
lar grain. Then, we compute the spinning dust emissivity 
and compare the results with those from disk-like grains. 

4.1. Power Spectrum 

A torque-free rotating irregular grain can rotate either 
in the positive or negative flip state, and the treatment in 
this section is applicable for a general flip state. To flnd 
the electric dipole emission spectrum from such an iso- 
lated irregular grain, we follow the approach in HDLIO. 
We first represent the dipole moment fj, in an inertial co- 
ordinate system, and then compute its second derivative. 
We obtain 

3 

where fii are components of /x along principal axes a,, 
ai are second derivative of a^ with respect to time, and 
i = 1,2 and 3. 

The instantaneous emission power of the dipole mo- 
ment is defined by 

PMq,t) = ^ii''. (10) 

The power spectrum is then obtained from the Fourier 
transform (FT) for the components of jl. For example, 
the amplitude of fix at the frequency Vk is defined as 

/ + 00 
jlx{t)exp {-i2Trv'kt) dt, (11) 
-oo 

where k denotes the frequency mode. The emission 
power at the positive frequency Uk is given by 

fed,fe( J, q) = A(/i^,fe + ij-lk + i^lk), (12) 



where the factor 2 arises from the positive/negative fre- 
quency symmetry of the Fourier spectrum. To reduce the 
spectral leakage in the FT, we convolve the time func- 
tion fj, with a Blackman-Harris window function (Harris 
1978). The power spectrum needs to be corrected for the 
power loss due to the window function. 

The total emission power from all frequency modes for 
a given J and q is then 




where T is the integration time. ^ 

Here we compute power spectra for case 1 (/xi = /x/ a/S) 
of fj, orientation, but our approach can be used straight- 
forward for any orientation of /u. 

Figure 2 presents normalized power spec- 
tra (squared amplitude of Fourier trans- 
forms), |FT(^,,j,)|Vmax(|FT(/i^)|2) and 
|FT(^2)P/max(|FT(/ix)P) for the components jl^ 
(or fly) and jl^, for a torque- free rotating irregu- 
lar grain having the ratio of moments of inertia 
Ii : I2 : I3 = 1 : 0.6 : 0.5 in the positive flip state with 
various q. Circles and triangles denote peaks of the 
power spectrum for components of fix (or fly) and flz, 
respectively. 

Multiple frequency modes are observed in the power 
spectra of the irregular grain, but in Figure 2 we show 
only the modes with power no less than 10^"^ the maxi- 
mum value. Both positive and negative flip states have 
the same frequency modes. 

For q < /1//2, we found that power spectra for fix (or 
fly) have angular frequency modes 

Wn^«(</')+m(|V'|), (14) 

where the bracket denotes the averaging value over time, 
and m = 0, ±1, ±2... denote the order of the mode. The 
frequency modes for jlz arc given by 

LOn = nm, (15) 

where n is integer and n > 1. 

The above frequency modes arc equivalent to results 
obtained from a quantum mechanical treatment. For the 
limiting case of an axisymmetric grain with Ii > I2 = I3, 
the grain rotational Hamiltonian is H = J^/2J2 — (/i — 
/2)/(2/i/2) where Jz = JcosO is the projection of J 
on the symmetry axis. From a quantitization, — > 
J{J + l)h^ and Jz Kh where J and K are quantum 
numbers for the angular momentum and its projection on 
the grain symmetry axis, the rotational energy is writ- 
ten as Ej^k = J{J + l)h^ /^h - K^h^h - h)l {^hh). 
For PAHs, J is usually very large, so an emission mode 
arising from the transition from the energy level Ej^k to 
Eji^K' have frequency given by 

1 fdEjK . dEjK , \ 

where BEj^k/OJ « {Jh^/h) = H and dEj.x/dK = 

{Kh'^/Iil2){Ii — h) = fi'tp- The selection rules for elec- 
tric dipole transitions correspond to AJ = 0, ± 1, 

^ This is the result of Parseval's Theorem. 
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and /S.K = 0, ± 1, so the total emission modes are 
ijJjj',KK' = <j>, (p i ip and ip, as found in HDLIO. The 
rotational energy levels of asymmetric rotors are speci- 
fied by three quantum numbers - the total angular mo- 
mentum J and two additional quantum numbers K^i 
and K+i. Electric dipole transitions obey selection rules 
A J = 0,±1, AK^i = ±1.±3, and AK+i = ±1.±3 
(see, e.g., Townes & Schawlow 1955) ^ . In the large J 
limit, these selection rules will correspond to the modes 
of Equations (14) and (15). 

For q > /1//2, only the mode to = has angular fre- 
quency given by Equation (14). 

In the following, the emission modes induced by the 
oscillation of Hx or /Xy, which lie in the xy plane, per- 
pendicular to J, are called in-plane modes, and those 
induced by the oscillation of /Xz in the direction perpen- 
dicular to the xy plane, arc called out-of-plane modes. 
The order of mode is denoted by m and n, respectively. 

Figure 2 also shows that the emission power for out- 
of-plane modes a;„ is negligible compared to the power 
emitted by in-plane modes w^. 

The upper panel in Figure 3 shows the variation of the 
angular frequency of the dominant modes to = 0, ±1, ±2 
with q. The kink at g = (jsp — I1/I2 in LOm arises from 
the change of grain rotation state. For q < qgp, the 

frequency oJm^o = {4) is almost constant, whereas it 
increases with q to Wm=o = J/H a-t q = Ii/I^. The 
frequency of the mode m = 1 also increases with q for 
q > qsp. The increase of frequency for the two dominant 
modes to = and 1 for g > q^p can result in an increase 
of peak frequency of spinning dust spectrum. This issue 
will be addressed in the next section. 

The middle panel in Figure 3 shows the emission power 
at the dominant frequency modes |to| < 2. For q ^ 1, 
most of emission comes from the mode m — —1. As q 
increases, the mode to = becomes dominant. 

The total emission power from all modes for irregular 
grains of different : Iz and 7i : 72 = 1 : 0.6 is shown in 
the right panel. 

We compute emission power spectra to infer Wfc and 
Ped.fe, as functions of J and g, for various ratio of mo- 
ments of inertia /i : /2 : -^3- The obtained data will be 
used later to compute spinning dust emissivity. 

4.2. Spinning Dust Emissivity 

To calculate the spinning dust emissivity, the first step 
is to find the distribution function for grain angular mo- 
mentum. 

4.2.1. Angular Momentum Distribution Function 

Following HDLIO, to find the distribution function for 
the grain angular momentum J, we solve Langevin equa- 
tions describing the evolution of J in time in an inertial 
coordinate system. They read 



dJi = Aidt -\- 



for i 



(17) 



X, y, z, 

where dqi are random Gaussian variables with {dq^) = 

dt, Ai = (AJi/At) and Bu = ((AJj)^ /At) are damping 
and diffusion coefficients defined in the inertial coordi- 
nate system. 

^ The rotational levels of a particulaj: asymmetric rotor, H2O, 
are shown in Draine (2011). 



For an irregular grain, to simplify calculations, we 
adopt the Ai and Bu for a disk-like grain obtained in 
HDLIO. Following DL98b and HDLIO, the disk-like grain 
has radius R and thickness L = 3.35 A, and the ratio of 

moments of inertia along and perpendicular to the grain 
symmetry axis h = Thereby, the effect of non- 

axisymmetry on Ai and Ba is ignored, and we examine 
only the effect the wobbling resulting from the irregular 
shape has on the spinning dust emissivity. 
In dimensionless units, J' = J//[|a;x.|| with ujt.\\ = 

1 /2 

(2fcBTgas//||) being the thermal angular velocity of 



the grain along the grain symmetry axis, and t' 
Equation (17) becomes 

dJl = A'.dt' + 

where {dq^^) = dt' and 

J' 2 



where 



7h = 

'''ed,eS 



^gas,eff ^ed.eff 



Bi, 



2/11 fcRT. 



-■^H,||, 



B-*^ gas 



tot, II 



th.II cos^ ^ -I- 7h sin^ 6* ' 

-Ptot,-LTH,|| 
-Ptot,||TH,_L ' 
Ted.eff 
^H,|| ' 



(18) 

(19) 
(20) 

(21) 

(22) 
(23) 



where Thji and th^j. are rotational damping times due to 
gas of purely hydrogen atom for rotation along parallel 
and perpendicular direction to the grain symmetry axis 
ai) Ted,e« is the effective damping time due to electric 
dipole emission (see Appendix B), ^ is the angle between 
ai and J, and -Ftotj and Ftot,_L are total damping coeffi- 
cients parallel and perpendicular to ai (sec HDLIO). 

At each time step, the angular momentum Jj obtained 
from LEs is recorded and later used to find the distribu- 
tion function fj with normalization fjdJ = 1. 

4.2.2. Emissivity from Irregular Grains 

As discussed in §4.1, an irregular grain rotating with a 
given angular momentum J radiates at frequency modes 
ujk = LOm with to = 0, ±1,±2... and uok ^ with n = 
1, 2, 3... (see Eqs 14 and 15). For simplicity, we define the 
former as and the latter as where i indicates the 
value for to and n. These frequency modes depend on 
the parameter q{s), which is determined by the internal 
thermal fluctuations within the grain. 

To find the emissivity at an observation frequency v = 
uj/2tt, we need to know how much emission each mode 
LUk contributes to the observation frequency ly. 

Consider an irregular grain rotating with the angular 
momentum J, the probability of finding the emission at 
the angular frequency w depends on the probability of 
finding the value oj such that 



pdf{uj\J)du; = /vre(s, J)ds = Aexp I - ^^^^^^ ^ 



ds, (24) 
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Fig. 2. — Normalized power spectrum of a torque-free rotating irregular grain with Ii : I2 ■ I3 = 1 : 0.6 : 0.5 for different values of 
q = 1.05,1.60 (i.e. q < qsp = h/h) and q = 1.81 > Qsp. The components of |FT{/i^)|V max(|FT(/i^)|2)) (or |FT(/ia)p/ max(|FT(/ia;)P)) 
are indicated by circles, while the components of |FT(/i2)p/ max(| FT (/ia; ) p ) are indicated by triangles. Orders of in-plane modes m and 
out-of-plane modes n are indicated, and case 1 (fii = /i/\/3) of orientation is assumed. 
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Fig. 3. — Left: Angular frequency of in-plane modes as a function of q for an irregular grain with Ii : I2 : I3 = 1 : 0.6 : 0.5. The kink at 
q = qsp = I1/I2 corresponds to the change from rotation around ai to rotation around 33. Middle: Emission power radiating at frequency 
modes ujm as a function of q. Right: Total emission power from all frequency modes as a function of q for various ratios I2 : I3 while 
II : I2 = ^ : 0.6 is kept constant. Case 1 [fii = tJ./V3) of /X orientation is considered. 



where we assumed the VRE regime with /vre given by 
Equation (8). 

For the mode lu = Wfe(s), from Equation (24) we can 
derive 



pdfk{uj\J) 



f dujk \ 



/vre(s, J). 



\ ds J 

The emissivity from the mode k is calculated as 
1 /•^" 

= ^ PcdAJ.<l<)fj{J)pdfk{co\J)27TdJ 



(25) 



1 

in 



PedAJ,q>)fj{J)pdfkHJ)2Tr dJ, (26) 



where q< and q> denote q < qsp and q > qsp, respec- 
tively; Ji and Ju are lower and upper limits for J corre- 
sponding to a given angular frequency u!k{J,q) = 1^, and 
27r appears due to the change of variable from i> to lo. 

Emissivity by a grain of size a at the observation fre- 
quency v arising from all emission modes is then 



' 7 ,Jv,k- 
k 



(27) 



Consider for example the emission mode k = mp. For 
the case I2 slightly larger than I3, this mode has the an- 
gular frequency iVmo = {4>) — {J/h)qo with qo obtained 
from calculation of ujmn, is independent on q for q < gsp-^ 

^ go approaches I1/I2 as I3 — > I2, i.e., when irregular shape 
becomes spheroid. 



As a result 



pdfmo{^\J) =S{uj- (J//i)go) ■ 



(28) 



Thus, the first term of Equation (26), denoted by j^^^ ^, 
is rewritten as 



1 



Pcd,„i„(J, q<)fj{J)S (w - {J/h)qo) dJ, 

Ji 

1 hfjjJ o) 

2 qo 



--Pcd,ino(>/O,g(s)), 



(29) 



where Jq — Iiuj/qo, and the value of q{s) remains to be 
determined. 

For q > qsp, {(p) is a function of q. Hence, the emissivity 
(26) for the mode k = irq becomes 
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dsPcd,m„iJo,q{s))fYHEiJo,s) 



Pod,mo(J, q>)fjiJ)pdUMJ)dJ, (30) 



where Sgp is the value of s corresponding to q = qsp, 
and the term Pcd.mo('A)j (/(s)) in Equation (29) has been 
replaced by its average value over the internal thermal 
distribution /vre- 

The emissivity per H is obtained by integrating over 
the grain size distribution: 



1^ = ± 
where is given by Equation (26) 



da 



(31) 



4.2.3. Emissivity from Disk-like Grains 

HDL 10 showed that a disk-like grain with an angular 
momentum J radiates at frequency modes 



^,=<P + iip = — [h + i{l-h)cose], (32) 



C0ni=1p = -r(l - h)cOS0, 
-(11 



(33) 



where z = and ±1. 

The emission power of these modes are given by fol- 
lowing analytical forms (HDLIO; Silsbee et al. 2011): 



6c3 

^< sm^, 



(34) 
(35) 

(36) 



For the disk-like grain, from Equation (5), we obtain 
the number of state in phase space s for q spanning from 
1-q: 



s = 1 



h — q 



1/2 



1 — cos 9, 



(37) 



h-1 

where q = 1 + (h — l) sin^ 9 has been used. Thus, for an 
arbitrary mode with frequency ojk, we obtain 

pdfk{ui\J)doj = /vre(s, J)ds = /vre(6', J) iim9d9. (38) 

Taking use of a; = a;fe( J, ^), we derive 

pdfkHJ) = /vre(^, J) sin^. (39) 

Therefore, by substituting Equations (34)- (36) in 
Equation (26), the emissivity at the observation fre- 
quency v = lj/{2'k) from a disk-like grain of size a is 
now given by 



3c3 



w (sm 



2 6c3 . 


/■'■ 


iJi 


+ 26?'" . 


/" 


iJi 


2 3c3 . 


r 


i.h 



pdfnMJ)fj{J)dJ, (40) 



where pd/ni±i and prf/m arc easily derived by using Equa- 
tion (39) for Wm±i and Wm, and J; = I^^uj/{2h— 1) and 
Ju = I\\uj for m±i mode, and Ji — I\\io/{h — 1) and 
Ju = oo for ni mode (see Eqs 32 and 33) . 

As in HDLIO, wc also consider the regime with- 
out internal relaxation. Then, the distribution func- 
tion /vre(^, J) is replaced by a Maxwellian distribution 
(hereafter Mw regime): 

/M^(g) = A 1 (41) 



^'^ {cos^9 + hsm^ 9) 



For symmetric rotors {h = 1), the distribution function 
for 9 is isotropic: 



h-o{9) - ^, 



(42) 
condition 



satisfying the normalization 

/o''/iso(^)27rsin^d^ = 1 

4.3. Results 

In this section, we present results for spinning dust 
emissivity from grains of irregular shape obtained using 
the method in Section 4.2. 

4.3.1. A Simple Model of Irregular Grain 

We assume that smallest grains of size o < 02 = 6 A, 
have planar shape, and larger grains are spherical. To 
compare with emission from the disk-like grains, we con- 
sider the simplest case of the irregular shape in which the 
circular cross-section of the disk-like grain is changed to 
elliptical cross-section. We are interested in the emission 
of two PAHs with the same mass M and thickness L, 
therefore, the length of the semi-axes of elliptical disk is 
constrained by 



M 



(43) 



where R = (4a3/3L)^/^ is the radius of the disk-like 
grain, 62 and 63 arc the length of semi-axes a2 and 33, 
and fei = L is kept constant. Assuming that the circu- 
lar disk is compressed by a factor a < 1 along a2, then 
Equation (43) yields 



aR, 63 



(44) 



Denote the parameter 77 = 63/62 = a , then the degree 
of grain shape irregularity is completely characterized by 

The moments of inertia along the principal axes for the 
irregular grain are 



(45) 



where is the moment of inertia for the disk-like grain. 
The ratios of moment of inertia read 



h 
h 



3(r? 



■n 



h 
h 



{L/Rf 
{L/Rf- 



^377' 
3r/ 



(L/i?)2 + 3?7-i' 



(46) 
(47) 



For each grain size a, the parameter rj is increased from 
77 = 1 to ry = T/max- Howevcr, rymax is constrained by the 
fact that the shortest axis a.2 should not be shorter than 
the grain thickness L. We below assume conservatively 

^max ^ 3/2. 

4.3.2. Rotational Emission Spectrum 

Although irregular grains radiate in large number of 
frequency modes, we take into account only the domi- 
nant modes with order \m\ < 2, and disregard the modes 
that contribute less than ^ 0.5% to the total emission. 
We also assume that grains smaller than 02 have a fixed 
vibrational temperature Tvib, and that for the instanta- 
neous value of J the rotational energy has a probability 
distribution /vre (i-e. VRE regime). 
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Wc adopt the grain size distribution dn/da from Draine 
& Li (2007) witli the total to selective extinction i?v = 
3.1 and the total carbon abundance per hydrogen nucleus 
be = 5.5 X 10~^ in carbonaceous grains with amin = 
3.55 A and flmax = 100 A. 

The spinning dust emissivity is calculated for a so- 
called model A (similar to DL98b; HDLIO), in which 
25% of grains have the electric dipole moment parame- 
ter /3 2/3o, 50% have /3 = /3o and 25% have /3 = 0.5/3o 
with (3o = 0.4 D. In the rest of the paper, the notation 
model A is omitted, unless stated otherwise. 

Figure 4 shows the spinning dust spectra for different 
degree of irregularity rj and various dust temperature Tvib 
in the WIM. The emission spectrum for a given Tvib shifts 
to higher frequency as r] decreases (i.e. the degree of grain 
irregularity increases), but their spectral profiles remain 
similar. 

To see more clearly the effect of grain shape irregular- 
ity, in Figure 5, we show the variation of peak frequency 
i^peak and peak emissivity jpeak as a function of r] for dif- 
ferent Tvib and for case 1 (/Ui = n/V^) of /2 orientation 
for the WIM. 

For low Tvib (i-e. Tvib = 20 K), fpeak increases slowly 
with T], while jpeak varies slightly. As Tvib increases, both 
fpeak and jpeak increase rapidly with r], but the former 
increases faster. 

Figure 6 shows the results for case 2 {^i = 0) of /x 
orientation. Similar to case 1, t'poak increases with rj for 
Tvib > 60 K, but changes rather little for Tvib = 20 K. 
The variation of jpeak with rj is more complicated. For 
Tvih < 60 K, it increases first to ~ 1.23, and then drops 
as T] increases. 

5. SPINNING DUST SPECTRUM: EXPLORING 
PARAMETER SPACE 

Below we explore parameter space by varying dust 
temperature, the lower grain size cutoff, the grain dipole 
moment, and the gas density. Our results from previous 
section indicate that the spinning dust spectrum from 
grains of irregular shape is shifted in general to higher 
frequency, relative to that from disk-like grains, but their 
spectral profile remains similar. Therefore, it is sufficient 
to explore parameter space of spinning dust for disk-like 
grains, extrapolating to irregular grains using our results 
from the previous section. Spinning dust emissivity is 
then calculated using Equation (31) in which j° is given 
by Equation (40). 

5.1. Effect of internal thermal fluctuations 

In this section, we account for the fluctuations of the 
dust temperature in smallest grains. 

Let Td be the decoupling temperature for the V-R en- 
ergy exchange. For Tvib > 2d, the V-R energy exchange 
is efficient, and Trot ~ ^vib- For the instantaneous value 
of J the rotational energy has a probability distribu- 
tion (8) determined by V-R energy exchange (the VRE 
regime). For Tvib < 7d, no V-R energy exchange exists. 
The rotational emissivity averaged over temperature fluc- 
tuations for a fixed grain size reads 

/dP 
j^(Tvib > Ta)^dT,i^ 

+ I fATrot < Ti)PoS{T,ot - Td)rfTrot, (48) 



where 

Po = / ;^rfTvib, (49) 

Jo "-tvib 

is the probability of the grain having Tvib < ?d, and 
dP/dT is the probability for the grain having tempera- 
ture mT,T + dT. 

Figure 7 shows spinning dnst spectra averaged over 
temperature distribution, {jv{T^))T , for various decou- 
pling temperature Td compared to spectra for grains hav- 
ing a single vibrational temperature Tvib = denoted 
by ji/(Tvib = Td). For low Td < 70 K, the averaged emis- 
sivity {jv{Td))T is larger than ju{T^ih = Td), but the dif- 
ference between j,y(Tvib = Td) and {jv{Td))T decreases 
as Td increases. 

For Td = 60 K, Figure 7 shows that {j^{Td))T is ap- 
proximate to j,y(Tvib = Td) with the diff'erence less than 
a few percent. In the rest of the paper, we adopt a con- 
servative value Td = 60 K and calculate the spinning 
dust emissivity assuming that all grains have a single 
temperature Td. 

5.2. Minimum size ami„ 

The spinning dust emission spectrum is sensitive to 
the population of smallest dust grains, and its peak 
frequency is mostly determined by the smallest grains. 
When amin is increased, the peak frequency fpeak de- 
creases accordingly. 

Figure 8 shows the variation of fpeak as a function of 
flmin for various environments, and case 2 (/./i = 0) of /x 
orientation and VRE regime (Td = 60 K). As expected, 
fpeak decreases generically with amin increasing. 

5.3. Effect of electric dipole moment 

5.3.1. Variation of orientation of electric dipole moment 

Silsbce ct al. (2011) found that the total emission 
power Ped and the peak frequency are rather insensitive 
to the orientation of the electric dipole moment. That 
stems from their calculations assuming no internal relax- 
ation (i.e. Tvib oo), so that the angle 6 is drawn from 
a isotropic distribution function fiso- Below we study the 
variation of emissivity with the orientation of dipole 
accounting for the internal relaxation. 

To investigate the variation of ji, with the orientation 
of /IX, in addition to case 1 (/Ui = /i/\/3) and case 2 [fjbi = 
0) of /X orientation, we consider values oi y?^/ between 
2/3 (case 1) and 1 (case 2). We also consider the regime 
of fast internal relaxation (VRE) and without internal 
relaxation (Mw and iso). Model A for the distribution 
of electric dipole moment is adopted. 

Figure 9 shows the variation of t'peak [upper panel) 
and jpeak {lower panel) of spinning dust spectrum with 
It can be seen that for the Mw and iso f{9), t'peak 
changes rather slowly with fi'^, with somewhat larger 
variation found in the VRE regime. The peak emissivity 
for all regimes exhibits similar trend that decreases with 
increasing n"^, though such a decrease is only ~ 10% 

when n orientation changes from case 1 {jjLi = fJ,/V^) to 
case 2 (/xi = 0). 

5.3.2. Variation of distribution of electric dipole moment 
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Fig. 4. — Emissivity per H from irregular grains of different degree of irregularity r\ = 63/^2 with Tvib = 20, 60 and 200 K for case 1 
{fjL-i = /x/\/3) of dipole moment orientation in the WIM. The emission spectrum shifts to higher frequency as r] decreases. The grain mass 
is held fixed as r] changes. 
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Fig. 7. — Emissivity per H for case 2 (^i = 0) of the /i orienta- 
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ity spectra averaged over temperature distribution, {ju{T^))T) a^nd 
spectra for grains having Tvib = 7d, jjy(Tvib = Tj), respectively. 
Two values = 20 and 60 K are chosen for illustration. 



Fig. 9. — Variation of the peak frequency t'peak (upper) and peak 
emissivity (lower) with orientation ( (/Xj + Ma ) /a*^ ) for VRE, Mw 
and iso regimes. Vpeak for the VRE regime decreases faster than 
that for the Mw and iso regimes as fJ, orientation changes from 
case 1 (/ii = ii/\/3) to case 2 (fj,i = 0), but the decrease is within 
10%. 
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Fig. 8. — Decrease of the peak frequency J^peak of spinning dust 
spectrum with the lower cutoff of the grain size distribution amin 
for various environment conditions. 



Our standard model (model A, see Table 1) assumed 
grains of a given size a to have three values of dipole mo- 
ment, (3 = /?o, /3o/2 and 2/3o with (3q = 0.4 D, appearing 
in Equation (2). 

Here we examine the possibility of considerable vari- 
ation in the dipole moment per grain, by considering a 
model (model B) in which 40% of grains have dipole mo- 
ment parameter /3 = /3o; 15% have /? = 2/3o and /3o/2; 
10% have /3 = 4/3o and ;3o/4; and 5% have /3 = 8^o and 
/3o/8. Emissivity per H for this model is compared with 
the typical model A in Figure 10. 

Figure 10 shows emission spectra for model A and B in 
cases 1 (/ii = /i/V3) and 2 (/ii = 0) of /i. orientation, for 
the CNM and WIM. Resuhs for the VRE (Td = 60 K) 
regime are presented. The peak frequency from model 
B is lower than for model A, but model B has a higher 
^■^ peak emissivity. The larger value of {jP) = 3.12/3q in 
model B leads to stronger electric dipole emission but 
also stronger rotational damping. 

5.4. Variation of values of dipole moment and gas 
density 
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Fig. 10. — Emissivity per H for the CNM (upper) and WIM 
(lower) for model A and B of dipole moment distribution and for 
case 1 (^1 = fi/\/3) and case 2 (^i = 0) of /X orientation and for 
the VRE (Td = 60 K) regime. 

Here we assume that all physical parameters are con- 
stant, including gas temperature and dust temperature. 
Only the value of characteristic dipole moment /3o and 
gas density nn are varied in calculations. We run LE 
simulations for 16 values of (3o from 0.05 — 3.2 D, and 
32 values of nn from 10^^ — Icm^'^ for the WIM and 
from 1 — 100 cm~^ for the CNM. For a given value of 
/3o, the gaseous rotational damping time and the damp- 
ing and excitation coefficient from infrared emission vary 
with the gas density nn as 



(50) 



Fm^F,n{nK) — , Gir = GiR(nH) — , (51) 

where tTh is the typical gas density given in Table 2. 
Other rotational damping and excitation coefficients are 
independent of gas density. 

The obtained distribution function for grain angular 
momentum is used to calculate spinning dust emissivity 
as functions of /3o and nu (see Sec. 4.2.3). 

Figure 11 shows the contour of peak frequency I'pcak 
and peak emissivity jpcak in the plane of njj, /3o for case 2 
(yUi = 0) oi fi orientation in the CNM and WIM and for 



the VRE regime with Td = 60 K. For a given /3o, both 
i^pcak and jpoak incrcasc with uh due to the increase of 
coUisional excitation with nn- For a given nu, as /3o 
increases, t^pcak decreases, but jpeak increases. That is 
because as /3q increases, electric dipole damping rate in- 
creases, which results in lower rotational rate of grains. 
Meanwhile, the emissivity increases with /3o as /3q. 

The results (i^pcak and jpeak) for the Mw regime are 
slightly larger than those for the former one due to the 
lack of internal relaxation. 

From Figure 11 it can be seen that the amplitude of 
variation of i^pcak and jpoak as functions of nn and /3o is 
very large. Thus, for a fixed grain size distribution, nn 
and /3o are the most important parameters in character- 
izing the spinning dust emission. 



6. 



EFFECT OF TURBULENCE ON SPINNING DUST 
SPECTRUM 



Turbulence is present at all scales of astrophysical sys- 
tems (see Armstrong et al. 1995; Chepurnov & Lazarian 
2010). The effects of turbulence on many astrophysi- 
cal phenomena, e.g. star formation, cosmic rays trans- 
port and acceleration, have been widely studied. In the 
context of electric dipole emission from spinning dust, 
compressible turbulence will produce density variations 
that can affect grain charging, rotational excitation, and 
damping. We expect enhancement of spinning dust emis- 
sivity per particle from denser clumps, with an increase 
in the emission integrated along a line of sight. This 
issue is quantified in the following. For simplicity, we 
disregard the effect of turbulence on the grain charging. 

6.1. Numerical Simulations of MHD turbulence 

Gas density fluctuations in a turbulent medium is 
generated from Magneto-hydromagnetic (MHD) simula- 
tions. The sonic and Alfvenic Mach number are defined 
as usual Mg = (^), Ma — {^) with 6V are amplitude 
for injection velocity, Cs and va = Bq/^/Attp are sound 
speed and Alfvenic speed. 

The simulations are performed by solving the set of 
non-ideal MHD equations, in conservative form: 

|^-KV.(pv)=0, (52) 



V • 



Pvv+lP+|^)l~^BB =pf, (53) 



— - V X (v X B) = 0, (54) 

with V • B = 0, where p, v and p are the plasma density, 
velocity and pressure, respectively, I is the diagonal unit 
tensor, B is the magnetic field, and f represents the ex- 
ternal acceleration source, responsible for the turbulence 
injection. Under this assumption, the set of equations is 
closed by an isothermal equation of state p — cj.p. The 
equations are solved using a second-order-accurate and 
non-oscillatory MHD code described in detail in Cho & 
Lazarian (2003). Parameters for the models of MHD sim- 
ulations are shown in Table 3, including the sonic Mach 
number Ms and Alfvenic Mach number Ma- 
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Fig. 11. — Contour of peak frequency and peak emissivity in the plane of rtjj and /3o for the CNM (upper) and the WIM (lower), for 
the VRE (Td = 60 K) regime. Case 2 (fii = 0) of /X orientation is considered. Peak frequency increases with njj but decreases with /3o. 
Emissivity increases with both nj] and /3o. 



TABLE 3 
Models of 512^ MHD simulations 



Model 




Ma 


Va/Cs 




1 


2.0 


0.7 


2.9 


1.7 


2 


4.4 


0.7 


6.3 


2.1 


3 


7.0 


0.7 


10 


3.0 



E 
o 



10.0 



1.0 



6.2. Results: influence of Gas Density Fluctuations 

In a medium with density fluctuations, the effective 
emissivity is 

r-l 



0.1 



0.) 



(55) 



model A 
case 2 
VRE, T, 



where f{x)dx is the fraction of the mass with p/{p) G 
{x,x + dx). We use compression distributions f{x) ob- 
tained from MHD simulations for Ms = 2 and 7 to evalu- 
ate for the WIM and CNM, respectively. We assume 
case 2 (/ii — 0) oi fi orientation. It can be seen that the 
turbulent compression increases the emissivity, and shifts 
the peak to higher J^pcak- 

The increase of peak frequency and peak emission in- 
tensity, as a function of Ms are shown in Figure 13. 
(jpeak) and jpeak are emissivity for the turbulent and 




V (GHz 



Fig. 12. — Spinning dust emissivity per H in the presence of 
compressible turbulence with sonic Mach number Ms = 2 and 7, 
compared to that from constant gas density njj = "h for the CNM 
(red) and WIM (green). The peak emissivity is increased, and 
the spectrum is shifted to higher frequency due to compressible 
turbulence. Case 2 (fii = 0) oi fJ, orientation is considered 



uniform media, respectively, (t^pcak) 
frequency for the former and latter. 
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Fig. 13. — Increase of peak emissivity and peak frequency due 
to compressible turbulence as functions of Mach sonic number Ms . 
Results for both fast internal relaxation (VRE, = 60 K) and 
without internal relaxation (Mw) are shown. Case 2 (/ii = 0) of /Lt 
orientation is considered. 

and 7 are obtained from the spinning dust spectra for 
the WIM and CNM, respectively. For Ms = 4.4, results 
for the CNM are chosen. It is shown that the emission 
intensity can be increased by factors from ~ 1.2 — 1.4, 
and peak frequency is increased by factors from 1.1 — 1.2 
as Ms increases from 2 — 7. 

6.3. Observational constraining turbulence effects 

The discussion of interstellar conditions adopted in 
DL98 was limited by idealized interstellar phases. It is 
now recognized that turbulence plays an important role 
in shaping the interstellar medium. 

The distribution of phases, for instance, CNM and 
WNM of the ISM at high latitudes can be obtained from 
absorption lines. Similarly, studying fluctuations of emis- 
sion it is possible to constrain parameters of turbulence. 
In an idealized case of a single phase medium with fluc- 
tuations of density with a given characteristic size one 
can estimate the value of the 3D fluctuation by studying 
the 2D fluctuations of column density. A more sophisti- 
cated techniques of obtaining sonic Mach numbers^ have 
been developed recently (see Kowal et al. 2008: Esquivcl 
& Lazarian 2010; Burkhard et al. 2009, 2010). In par- 
ticular, Burkhart et al. (2010), using just column den- 
sity fluctuations of the SMC, obtained a distribution of 
Mach numbers corresponding to the independent mea- 
surements obtained using Doppler shifts and absorption 
data. With such an input, it is feasible to quantify the ef- 
fect of turbulence in actual observational studies of spin- 
ning dust emission. 

7. GALACTIC CMB FOREGROUND COMPONENTS AND 
FITTING MODEL 

7.1. Galactic CMB foreground components 

The Galactic CMB foreground consists of four com- 
ponents: synchrotron emission (J^''"), free-free emission 
(7*), thermal dust emission {II'^), and spinning dust 

* It may be seen that Alfven Mach numbers have subdominant 
effect on the distribution of densities (see Kowal, Laaaxian & Beres- 
nyaJj 2008). Thus in our study (see Table 3) we did not vary the 
Alfven Mach number. 



jsyn 



jtd 



jsd 



(56) 



The synchrotron emission consists of soft and hard 
components. The soft component has the spectrum 
/sy" ~ zy-/3 with p - 0.7 - 1.2, and is dominant be- 
low 1 GHz. The soft synchrotron is produced from elec- 
trons accelerated by supernova shocks that spiral about 
the Galactic magnetic field (Davies, Watson & Gutier- 
rez 1996). The hard component, also named "Haze", 
is present around the Galactic center, and its origin is 
still unknown (Finkbeiner 2004). The free- free emission 
arising from ion-electron scattering has a flatter spectral 
index ^ with /3 w 0.1, and is an important CMB 
foreground in the range 10— 100 GHz. The third compo- 
nent due to thermal dust emission (vibrational emission 
of dust grains) dominates above 100 GHz. Finally, the 
spinning dust component dominates the range 20 — 60 
GHz. 

Dobler et al. (2009) subtracted the synchrotron con- 
tribution using the Haslam et al. (i982) 408 MHz 
map plus a model for the "haze" component. The 
remaining emission was decomposed into two compo- 
nents: (1) a "thermal-dust-correlated" spectrum corre- 
lated with the SFD dust map (basically, 100 /im) and 
(2) an Ha-correlated spectrum correlated with observed 
(reddening-corrected) Ha. Both components included an 
"anomalous emission" component, peaking around 30-40 
GHz, that is attributable to spinning dust. 

The spinning dust intensity is given by 



rsd 



dl, 



(57) 



where ji, is the spinning dust emissivity, and the integral 
is taken along a line of sight. Below, wc constrain the 
physical parameters of spinning dust by fitting the model 
of foreground emissions to Ha-correlated and thermal 
dust-correlated emission spectra. 

PAHs are likely to be irregular, but wc do not attempt 
to determine the degree of irregularity in the present pa- 
per. Instead, we calculate spinning dust spectra for disk- 
like grains, and account for the effect of the grain shape 
irregularity by scaling the emission to account for irreg- 
ularity. Consider a simple model of the grain irregular 
shape, with the ratio of semi-axes r/ = 63/62 = 1.2 (see 
§4.3). Then, the model peak frequency and peak emissiv- 
ity are scaled with those from disk-like grains as follows: 



, ,mod , , . , X t'peak(r?) 

I'peak = i'peak(dlSk) X — — , 

t'pcakl.dlSKj 

Jpeak(??) 



ipTat = ipeak(disk)x 



(58) 

Jpeak(disk)' ^^^^ 



where we take t'peak(r/)/fpeak(disk) 
ipeak(??)/jpeak(disk) from Figures 5 and 6. 



and 



7.2. Fitting to Ha-correlated emission 

The Ha intensity is approximately given by (Draine 
2011) 



/hc«=0.361T4 



-0. 942-0. 0311nT4 



EM 



cm ° pc 



R, (60) 
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Fig. 14. — Fitting three-component model to Ho-correlated spec- 
trum: the contours in the plane of njj and /3o for VRE 
(T(j = 60 K) regime and case 2 (fii = 0) of fJ,. The parameters at 
which has minimum are shown. 
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Fig. 15. — A three-component fit model (solid line) to Ho- 
correlated spectrum using CMB5 data (diamonds) from Dobler et 
al. (2009) for VRE {upper) and Mw (lower) regimes using case 2 
(/il = 0) of fl. Three components of free-free emission, spinning 
dust emission and CMB bias are shown separately. As shown are 
best fits parameters and x^- 

where T4 — Te/(10'* K), Tg is the electron temperature, 
and EM is the emission measure. 
Dobler et al. (2009) have determined the spectrum 



jobs q£ ^j^g Ha-correlated emission in the WMAP five- 
year data. The thermal dust contribution has been al- 
most completely removed as part of the SFD- correlated 
component, and synchrotron emission has similarly been 
removed by its correlation with the Haslam et al. (1982) 
408 MHz map. Within the frequency range considered 
here, the residual thermal dust emission and synchrotron 
emission are assumed negligible. 

Thus, as in Dobler et al. (2009), we fit the Ha- 
correlated spectrum with three components: free-free 
emission, spinning dust emission, and CMB cross- 
correlation bias. We present the ratio of intensity to 
Ha intensity. The free-free emission is given by (Draine 
2011) 



Fn 



23 GHz 



kjy sr 



where 



i^O = 0.197rr"+°-9*21nT4^ 



(61) 



(62) 



The CMB bias is defined as dB^{T)/dT\T=2.7 k oc i^^ , 
and we obtain 



cmb 



-Co 



1 



23 GHz/ plc(i/) 



kJy sr 



(63) 



where Co is a free parameter, and pic is the Planck 
correction factor to convert thermodynamic temperature 
AT to antenna temperature at frequency v. 

The spinning dust component is taken as the emission 
from the WIM. We obtain 



rsd 



■ — Sdo X 



'Ha 



'Ha 



(64) 



WIM 



where [-f^'^Z-^Haj^jf^^ is calculated for our model for the 
WIM with a standard PAH abundance, and the fitting 
parameter Sdo allows for the possibility that the PAH 
abundance in the WIM may differ from our standard 
value. 

Using Equation (60) for the above equation, we obtain 
the spinning dust intensity per Ha 



sd 1 



'Ha 



WIM 



8.55 X lO^s 



T 



jHa / WIM 

cm^'^/nu 
(n-e/nn)^ 

0.942+0.031 lnT4 



X ■ 



R 



kJy sr" 



The model emission intensity per Ha is then 

mod 



^Ha 



= <Fr 



■»( 



-Sdo X 



23 GHz 



'Hct 



-0.12 



WIM 



c, 



1/ 



2 



1 



23 GHz/ plc(i^) 



(65) 



kJy sr 



(66) 



For the spinning dust model, we consider a range of 
gas density uh £ [0.01, 1] cm~'^ and the value of dipole 
moment /3o G [0.1, 1] D. For each (nn, /3o), we compute 
(ii^/?^H)wiM- then vary the values of Fq, 5* do and Co- 
From a given value Fo we derive the gas temperature 
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using Equation (62). When T4 and Oi//'T'h)wim 
known, we obtain [-^^'^/-fHaj^ijyi from Equation (65). 

The fitting process is performed by minimizing the 
function 



X 



I jn 

E- 



mod 



robs 1 2 



(67) 



where 1°^^ is the measured Ha-correlated emission in- 
tensity, and is the mean noise per observation at the 
frequency band ly (sec Hinshaw et al. 2007; Dobler et al. 
2009). 

In general, is a function of three amphtude parame- 
ters Fo, Co, Sdo, and two physical parameters nu and /3o- 
Therefore, the fitting proceeds by minimizing x^ over five 

parameters. 

Figure 14 shows x^ contours as functions of /3o and nu 
for case 2 = 0) and for VRE {T^ = 60 K), for the 
values of Sdo and Fq for which x^ is minimized. We can 
see that the distribution of is locahzed and centered 
around the standard values in the plane nn — /3o- 

Figure 15 shows the fit of the three-component model 
to the Ha-correlated foreground spectrum for VRE and 
Mw regimes for case 2 of /x. Both regimes provide a good 
fit to the data. 

Table 4 shows best fit parameters and for different 
cases of fi orientation and different regimes of internal 
relaxation. We can sec that both orientations of fi (case 1 
and case 2) can produce a good fit with low x^, but case 2 
exhibits a relatively better fit. The best fit gas density is 
in the range 0.08—0.15 cm~^. In case 1, the best fit dipole 
moment for the VRE and Mw regimes is /3q ^ 0.76 and 
1.0 D. In case 2, the best fit requires /3q ^ 0.65 D and 0.84 
D for VRE and Mw respectively. Case 1 requires higher 
/3o to reproduce the observations than case 2 because the 
effect of grain shape irregularity is more important for 
case 1 . Similarly, the Mw regime requires higher /3o than 
VRE because the effect of irregiilarity is stronger for the 
former case. 

From Table 4, it can be seen that the best fit value 
Sdo ~ 0.04 — 0.06 is significantly lower than the value 
Sdo = 0.3 in Dobler et al. (2009). Our lower Sdo stems 
from the higher value of f3o required, from the increase 
of emissivity in the HDLIO model compared to the DL98 
model used in Dobler et al. (2009), and the further (mod- 
est) increase in emissivity for irregular grains. As a re- 
sult, a higher depletion of small PAHs is required to ob- 
tain a good fit. 

In addition, from best fit values Fq in Table 4 we de- 
rive the gas temperature ^ 2700 — 3500 K as found 
by Dobler et al. (2009). This is lower than the typical 
temperature T^. = 8000 K usually assumed for the WIM. 
Dong & Draine (2011) proposed a model of three com- 
ponents that can explain the low gas temperature in the 
WIM. 

7.3. Fitting to thermal dust- correlated emission 

In addition to the Ha-correlated emission, the fore- 
ground induces a dust-correlated emission spectrum with 
a usual thermal dust emission component falling from 94 
to 40 GHz, and another component rising from 40 to 23 
GHz (see Bennett et al. 2003). The latter component is 
consistent with the spinning dust emission (DL98ab; de 



TABLE 4 

Best pit parameters and for Ha-CORRELATED EMISSION 





/int(e) 


nil (cm-3) 


00 (D) 




Fo 


Co 




case 1 


VRE 


0.15 


0.76 


0.04 


0.09 


0.0012 


1.9 




Mw 


0.08 


1.0 


0.06 


0.1 


0.001 


4.5 


case 2 


VRE 


0.11 


0.65 


0.06 


0.09 


0.0012 


1.1 




Mw 


0.14 


0.84 


0.04 


0.09 


0.0012 


1.6 



Oliveira et al. 2004; HDLIO). Although this peak fre- 
quency around 22 GHz is consistent with prediction by 
the DL98 model, it is lower than the prediction by the 
improved model of HDLIO, using the same parameters 
for the CNM as in DL98b. 

As in Dobler et al. (2009), we fit the thermal dust- 
correlated emission spectrum with a three-component 
model including spinning dust, thermal dust and CMB 
bias.^ 

Dobler et al. (1999) determined the spectrum 
/j^/T94GHz(FDS) for the contribution to the inten- 
sity from the thermal-dust-correlated emission, where 
T94GHz(FDS) is the antenna temperature at 94 GHz pre- 
dicted for the FDS dust model (Finkbeiner, Davis & 
Schlegel 1999). 

Planck collaboration (2011c) showed that the thermal 
dust emission can be approximated by 



rtd 



where Tt 



td 



17.6 K and e{u) = 



)-25(A/250 Atm)-^-«cm-2H . 
We model the spinning dust contribution as 



lf = Sdo(^) . 



(68) 
0.92 X 



(69) 



Thus, we seek to fit the Dobler et al. (2009) spectrum 

by 



rmod 



T94GHz(FDS) 



(94^)1'^''^^^"'-^ 



(jV/^h)cnm 



e{iy)Bi,{Ttd)]i^=9i GHz 



(70) 



where A, Sdo and Co are adjustable parameters. For 
each trial rifj and /3o, we adjust A, Sdo and Co to mini- 
mize x^- We consider different possible values of nu and 
/3o. 

Best fit parameters and x^ for different situations are 
shown in Table 5. In case 1 (^1 = /i/VS), the best fit 
corresponds to njj = 6.0 — 6.5cm^'^, and /3o ^ 0.94 — 
0.97 D for VRE and Mw regimes. In case 2 {/ii = 0), 
the best fit corresponds to nu = 7.0 — 7.5 cm~^ and ^0 ~ 
0.95 D for VRE and Mw regimes. The significance of 
fitting is low in both cases 1 and 2. Case 2 fits slightly 
better (x^ = 30 - 35) than case 1 (x^ = 54 - 56). 

Figure 16 shows the best-fit to the thermal dust- 
correlated spectrum obtained using case 2 models, for 
both VRE and Mw regimes. 

^ Free-free emission is not important for the thermal dust- 
correlated emission spectrum. 



17 



100 



3 



1^ 



- 10 



VRE, T^=60 K, case 2 
Sclo=0.91 



20 




40 

V (GHz) 



100 



Mw, case 2 



100 



3 




40 

V (GHz) 



100 



Fig. 16. — A three-component fit model (solid line) to ther- 
mal dust-correlated emission spectrum (diamonds) for VRE (Tj = 
60 K, upper) and Mw (lower) regimes using case 2 (fii = 0) of fl. 
Three components of spinning dust emission, thermal dust emis- 
sion and CMB bias are shown separately but solid, dashed and 
dashed-dotted lines. As shown is the best fit value of Sdo. 

We see that the fitting to thermal dust-correlated emis- 
sion has large ■ The reason is that the curvature of the 
model spectrum is larger than that of the observed spec- 
trum for frequencies 22 — 44 GHz. 



TABLE 5 

Best fit parameters and for thermal dust-correlated 
emission 





/int(e) 


"H(cm-^) 


/3o(D) 


Sdo 


A 


Co 




case 1 


VRE 


6.5 


0.97 


0.74 


184.3 


4.5 


54 




Mw 


6.0 


0.94 


0.63 


176.1 


4.7 


56 


case 2 


VRE 


7.5 


0.95 


0.91 


192.4 


2.9 


30 




Mw 


7.0 


0.96 


0.88 


184.2 


3.5 


35 



7.4. Classical vs. quantum treatment 

Ysard & Verstraete (2010, hereafter YVIO) questioned 
the validity of classical mechanics in the DL98 model for 
spinning dust emission. YVIO put forward a quantum- 
mechanical formalism for treating the rotational emis- 
sion. However, they only calculated spinning dust emis- 
sion for the rotational transition AJ = — 1 and AK = 0, 



which is induced by the oscillation of the dipole compo- 
nent along the grain symmetry axis. 

In our view, a quantum description of spinning dust 
is unnecessary, because, as shown in DL98b for spherical 
grain, the angular quantum number 



J 



IlUJ 



72 ('^ 
' 20 



5/6 



VlOO K 



1/2 



(71) 



which shows that even smallest PAHs have angular quan- 
tum number J ^ 1, therefore, the classical treatment 
should be valid. 

Miville-Deschenes et al. (2008) made an attempt to 
separate the spinning dust component from the Galactic 
foreground components (mostly synchrotron and anoma- 
lous emission) using both the WMAP intensity and po- 
larization data. The inferred spinning dust spectrum is 
presented in Figure 17 (diamond symbols). 

Ysard, Miville-Deschenes & Verstraete (2010) pre- 
sented a fit to observation data for spinning dust, which 
is extracted from WMAP data for regions with lati- 
tude h = 22.4° (the spinning dust spectrum extracted 
from WMAP data using the quantum mechanical ap- 
proach from YVIO). They assumed a lower cut-off of the 
grain size corresponding to the number of carbon atom 
Ar^i„ = 24 and 48 for the CNM and WNM, and dipole 
moment (3q — 0.3 D. Converting to the grain size using 
the usual relationship a — 10(iVc/410)-^/^ A, we obtain 
Omin ^ 3.88 and 4.89 A, respectively. However, since 
Oniin chosen for the WNM is much larger than that for 
the CNM, the peak frequency of its emission spectrum is 
expected to be much lower than that for the CNM, be- 
cause the peak frequency decreases with Omin increasing 
(see §5.2). This disagrees with results shown in Figure 5 
of Ysard et al. (2010), where the peak frequency of the 
WNM (t^pcak ^ 45 GHz) is much larger than that for the 
CNM (z^peak - 22 GHz). 

To see how well the improved model of spinning dust 
in HDL 10 fits to observation data for selected regions in 
Ysard ct al. (2010), we assume the same spinning dust 
parameters as in Ysard et al. (i.e. gas density, dipole 
moment I3q and Omin)- 

The total emission intensity from the CNM and WNM 
with column density 7Vh(CNM) and iVH(WNM) is given 

by 



T-mod 



iVH(CNM) 



■sd 



iVH(WNM) 



CNM 



■sd 



WNM 



where [jX /^t^hJ is the spinning dust emissivity per H 
calculated for the CNM and WNM, respectively (see 
Table 2). The fitting proceeds by minimizing (Eq. 
67) with iVH(CNM) and iVH(WNM) as free parameters. 
We found that the best fit is achieved at iVH(CNM) = 



3.5 X 10^" cm-^ and iVH(WNM) = 1.4 x 10^ 



The 



best fit is shown in Figure 17 and is also shown. We 
can see that our improved model can reproduce the ob- 
servational data for high latitude regions. 

8. DISCUSSION 

8.1. Relative importance of effects 

The present paper extends further the improved model 
of spinning dust emission in HDL 10 by including the 
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Fig. 17. — Comparison of our study with the anomalous spec- 
trum extracted from the observation data in YVIO. The sohd hne 
is the total cmissivity resulting from the CNM and WNM with 
AfH(CNM) = 3.5 X 1020 cm-2 and A^hCWNM) = 1.4 x lO^i cm-2. 
Squared symbols show observation data. All grains have the char- 
acteristic electric dipole moment /3 = /So = 0.3 D. Case 2 {m = 0) 
of is considered. 

following effects: irregular grain shape, dust tempera- 
ture fluctuations due to transient heating, orientation of 
dipole moment, and fluctuations of gas density due to 
compressible turbulence. 

8.1.1. Power Spectrum of Grains of Irregular Shape 

First of all, we find that a torque-free rotating irregular 
grain radiates at multiple frequency modes with angu- 
lar frequency and a;„. depending on a dimensionless 
parameter q, the ratio of the rotational energy to the 
rotational energy along the axis of major inertia. 

In case 1 (/iti = At/\/3) of /lx orientation, the mode 
ojm-i dominates the dipole emission for q close to 1 (when 
ai makes a small angle with J). As q increases, the 
mode ujmo increases accordingly, and become dominant. 
Physically, the mode uim=o arises from the precession of 
the dipole moment component parallel to ai axis about 
the angular momentum. Its emission power increases 
with the angle 6 between ai and J, and is maximum 
when ai is perpendicular to J, i.e., 6 = n/2 or q = qsp. 

In case 2 (/xi = 0, i.e. no dipole component along ai), 
the mode m = does not exist when the grain rotates 
about ai for q < Qsp- For q > qsp, the irregular grain ro- 
tates about as, therefore, the dipole moment component 
along as induces the frequency mode m = 0. 

In both case 1 and case 2, the modes m = ±1, which 
rcstilt from the oscillation of dipole moments along a2 
and as, always exist, and have power increasing with q 
up to g = qsp, then decrease in importance as q increases 
further. 

8.1.2. Increase of spinning dust emission with the degree of 

grain shape irregularity 

For both case 1 (/xi = ^/-v/S) and 2 (^i = 0) of /2 
orientation, we find that the peak frequency increases 
rapidly with increasing degree of grain shape irregularity 
77 (see §4.3). The underlying reason is that the deforma- 
tion of the grain from the disk-like shape allows the grain 
to rotate along the axis of minimum moment of inertia 



as, in addition to rotation about ai and a.2 axes. For a 
given angular momentum, the angular velocity along as 
is largest, so that the emission frequency of the irregular 
grain is in general higher than that of the disk-like grain 
with the same mass. 

Peak emissivity also increases with the grain shape ir- 
regularity for case 1. But for case 2, j/pcak decreases with 
T] for Tvib < 60 K, and starts to increase with r] as Tyib 
increases. 

When the internal thermal fluctuations are rather weak 
{Tj > Tvib) SO that the grain tends to rotate with the 
minimum rotational energy, 9 and q deviate slightly from 
9 = and q = 1, and the rotational emission spectrum 
of irregular grains differs only slightly from that of disk- 
like grains. In contrast, for strong thermal fluctuations 
{Tj 4C Tvib), the thermal fluctuations enable grains to 
spend an important fraction of time rotating with large 
q. Because the total power emission increases with q, 
the higher probability of rotating with larger q results 
in larger emissivity and peak frequency than disk-like 
grains. 

8.1.3. Orientation of dipole moment 

We investigated the effects of changing the orientation 
of the electric dipole moment /j,. We found (see Fig. 9) 
that as the orientation is varied from having a compo- 
nent parallel to ai (case 1, ^ii = /x/V3) to being entirely 
perpendicular to ai (case 2, /ii = 0), the peak frequency 
i^poak and the peak emissivity both decrease by ^ 10%. 
Given the sensitivity of the spectrum to other variables 
(mass distribution, grain shape, gas density, gas tem- 
perature), it does not seem possible to infer the orienta- 
tion of /J, from observations of the spectrum. 

8.1.4. Density Fluctuation due to Compressible Turbulence 

We identifled interstellar turbulence as another factor 

that influences emission spectrum. The spinning dust 
emissivity is determined by two principal processes: col- 
lisional and radiative damping and excitation. Local 
compression due to turbulence increases the coUisional 
excitation, that results in increased spinning dust emis- 
sivity. We found that the emission is increased by a factor 
from 1.2 — 1.4 as the sonic Mach number increases from 
Ms = 2 - 7. 

8.2. Constraints from Spinning dust Emission 

Recent studies showed that the correspondence of the 
DL98 model to observations can be improved by adjust- 
ing the parameters of the model. For instance, the five- 
year ( WMAP) data showed a broad bump with frequency 
at ^ 40 GHz in the Ha-correlated emission (Dobler & 
Finkbeiner 2008; Dobler et al. 2009). Using Ha as a 
tracer of the WIM, they showed that this bump can be 
explained using the DL98 model by varying either dipole 
moment or gas density of the WIM. 

By fitting the 3 components of synchrotron, free-free 
emission and spinning dust emission to Ha-correlated 
spectrum, Dobler et al. (2009) found that the spin- 
ning dust model with /3o = 0.15 D and density njj = 
0.15 cm~"^ could reproduce the observed spectrum, with 
Vpeak = 37 GHz. Our best fit values from the improved 
model correspond to nn = 0.08 — 0.15 cm~^ and /3o in 
the range of 0.6 — 1.0 D. 
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In addition to the Ha-corrclatcd emission, WMAP 
data show a thermal dust- correlated spectrum that de- 
clines from 23 GHz to 40 GHz. This peak at low fre- 
quency is lower than the HDLIO model prediction using 
the standard parameters for the CNM. 

Our results show that the thermal dust-correlated data 
can be fitted with spinning dust emission from the CNM 
of density riH ^ 6— 8cm^^ and /3o ~ 0.95 — 1.0 D. Similar 
to Dobler et al. (2009), we also found that the best-fit 
model did not provide a very good fit (high x^). The rea- 
son for that is the rotational spectrum is steeper than the 
observation data in frequency range 20 — 50 GHz. Dobler 
et al. (2009) suggested that the superposition of spinning 
dust spectra from different ISM phases along a sight-line 
would produce a flatter spectral slope. With our results 
for irregular grains, one can see that by averaging the 
rotational spectrum over various degree of irregularity rj 
with the fraction of irregular grains decreasing with rj 
(see §4), the obtained spectrum becomes shallower, that 
can improve the fit. 

8.3. Range of applicability of the model of spinning dust 

emission 

The model of spinning dust emission has been used to 
interpret the anomalous microwave emission in the gen- 
eral ISM (e.g. Finkbeiner 2004; Dobler & Finkbeiner 
2008; Gold et al. 2009, 2011; Planck Collaboration 
2011c), in star forming regions in the nearby galaxy NGC 
6946 (Scaife et al. 2010; Murphy et al. ' 2010) and in 
the Persus and Ophiuchus clouds (Gassasus et al. 2008; 
Tibbs et al. 2010; Planck Collaboration 2011a). Early 
Planck results have been interpreted as showing an emis- 
sion excess from spinning dust in the Magellanic Clouds 
(Bot et al. 2010; Planck CoUaboration 2011b). 

This paper together with HDLIO presents a compre- 
hensive model of spinning dust, accounting for non-disk- 
like ("irregular") grains, new electric dipole emission 
modes from torque-free rotation of irregular grains, and 
investigating the effects of variation in grain dipole mo- 
ment, gas density, and different regimes of vibrational- 
rotational mode coupling. We believe that apart from 
being an important CMB foreground, the spinning dust 
spectrum can become an important diagnostic tool to 
constrain physical properties of astrophysical dust (e.g. 
size distribution, shape, electric dipole moment and gas 
density) in various environments. 



9. SUMMARY 

The model of spinning dust emission is further ex- 
tended by accounting for effects of irregular grain shape, 
fluctuations of dust temperature, and effects of ISM tur- 
bulence. We consider both regimes of fast internal relax- 
ation and without internal relaxation. Our main results 
are as follows. 

1. Small grains of irregular shape radiate in general 
at multiple harmonic frequency modes. The rotational 
emission shifts to higher frequency as the degree of grain 
shape irregularity increases, but the spectral profile re- 
mains similar. The effect of the grain shape irregularity is 
more important for higher dust temperature or stronger 
internal thermal fluctuations. Depending on the irregu- 
larity parameter rj, peak frequency and peak emissivity 
can be increased by a factor of up to 1.4, relative to disk- 
likes grains of the same mass (see Figs. 5 and 6) for case 
1 {m = n/Vs). 

2. Fluctuations of dust temperature T^ih also increase 
the rotational emissivity relative to the emissivity for 
grains of a steady low temperature. 

3. Fluctuations of gas density and gas pressure due to 
compressible turbulence enhance both emission and peak 
frequency of spinning dust spectrum compared to that 
in uniform media. An increase in emission by a factor 
from 1.2 — 1.4 is expected as the sonic Mach number Mg 
increases from 2 — 7. 

4. Spinning dust parameters (e.g., gas density and 
dipole moment) are constrained by fitting the improved 
model to WMAP cross-correlation foreground spectra, 
including Ha and thermal dust-correlated spectra. We 
find a reduced PAH abimdance in the WIM {Sdo « 0.05) 
with dipole moment parameter /3o ~ 0.7 D. For the 
thcrmal-dust-correlatcd emission, we find a normal PAH 
abundance {Sdo « 0.9) and Po « 0.95 D. 

5. Our improved model also provides a good fit to 
WMAP data for selected regions at high latitude (6 = 
22.4°) obtained by Ysard et al. (2010). 
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APPENDIX 

A. ELECTRIC DIPOLE EMISSION FROM AN IRREGULAR GRAIN 
Al. Torque-free motion 

The effective size a of an irregular grain with volume V is defined as the radius of a sphere with the same volume 
V, i.e., 

a = (3V/47r)'/^ (Al) 

In general, an irregular grain can be characterized by a triaxial ellipsoid with moments of inertia Ii , I2 and I3 around 
three principal axes ai, a2 and as, respectively. Define dimensionless parameters so that the moments of inertia are 
written as 

= ctjisph, (A2) 

where /gph = (87r/15)pa^ is the moment of inertia of the equivalent sphere of radius o, and p is the mass density of 
the grain. 
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For a torquc-frcc rotating grain, its angular momentum J is conserved , while the angular velocity a; nutates and 
wobbles with respect to J. We can term the wobbling associated with the irregularity in the grain shape irregular 
wobbling, to avoid confusion with thermal wobbling (also thermal fluctuations) induced by the Barnett relaxation 
(Purcell 1979) and nuclear relaxation (Lazarian & Drains 1999b). 

A detailed description of the torque-free motion for an asymmetric top in terms of Euler angles 6, cj) and tp (see Fig. 
1) can be found in classical textbooks (e.g. Landau & Lifshitz 1976; see also WD03), and a brief summary is given 
below. 

Consider an ellipsoid with three principal axes ai , a2 , as and moments of inertia Ii > I2 > Is- 
Let define a dimensionless quantity 



(J2-/3)(g-l) 

(7i-72)(l- W/i)' 



(A3) 



where q = 2IiE^^j^ / ,P is the ratio of total kinetic energy to the rotational energy along the axis of major inertia ai. 
(a) For q < I1/I2, < I, the solution of Euler equations is 



u^i = ±i (^) dn(r), 
where cn, sn and dn are hyperbolic trigonometric functions, and r is given by 



= tJ 



{h-l2)(l-hq/h) 



hhh 



1/2 



(A4) 

(A5) 
(A6) 

(A7) 



and the sign ± in lji and CJ3 are taken the same. We denote the rotation with + and — sign as positive and negative 
rotation state. For q < /1//2, the grain mostly rotate about the axis of major inertia ai, while it rotates about as for 
q > h/h- 

The rotation period around the axis of major inertia ai is 



where F is the elliptic integral defined by 

F{e,m) = I d0{l-msm^e) 
Jo 

[b) For q > /1//2, angular velocities are given by 

-3=±i(W)'''dn(r), 



1/2 



where 



Rotation period for this case is given by 



T = tJ 



il2-h){q-l) 



-,1/2 



/1/2/3 



Pr = 4F(7r/2, k- 



(A8) 
(A9) 

(AlO) 
(All) 
(A12) 

(A13) 
(A14) 



(c) For q ^ I1/I3, Equation (A3) shows that fc^ — > 00, the rotation of the grain is about the axis near as. Thus, 
0)3 « J/I3, and wi ~ a;2 ~ 0. From Euler equations, we obtain 



Iidu) 



dt 

/2dw2 



- = W2W3 {I2 - h) ■ 



dt 



= wswi {I3 - h). 



(A15) 
(A16) 
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Substituting 0J3 = flo = J/h, the equations are rewritten as 



^=JloJ2(l-^), (A17) 



dt V I2 

^ = -fioJifl-T^), (A18) 



dt V ii. 

where Jk = Ik<^k- Taking the first derivative of Equation (A17), and using (A18), we obtain solutions for Ji and .72- 

Ji=Acosojt, (A19) 

1 /2 

j^ = -(l-^] sincot = — (1 - ^ ) sinwi, (A20) 

where = ^1 — ^1 — , and ^ is a constant of integration. Denote A/ ^1 — ^ j = aJ with a is a small 
parameter, then 

Ji=aJ\l-j-j coscot, (A21) 
r \ 1/2 



J2 = aJ(l--^) sin Lot. (A22) 



The value of a is found by using the relation 

9=^. (A23) 

where we have use the fact that at t = 0, J2 = 0. Substituting Ji and J3 in Equation (A24) and plugging it into 
(A23), we obtain 



Hence, 

- h/h 



(A26) 



As g = /1//3, then a = 0, i.e., Ji = J2 = and J3 = J. 

When the angular velocity components are known, we can infer the orientation of the grain axes in the inertial 
coordinate system using Euler angles: 

cos^=:^i^, tan^=:^. (A27) 
A2. Flip states 

For a given J, there are two sets of solution (± sign) for of the Euler motion equations (see Eqs A4-A6). It can be 
seen that for q < q^p, two rotation states ± correspond to > and uji < 0, i.e., ai.J > and ai.J < 0. We define 
these rotation states as positive flip state and negative flip state (also WD03; Hoang & Lazarian 2008). For q > q^p, 
the similar situation occurs with ^3, and there are positive and negative flip states with respect to a3. 

A3. Electric Dipole Emission for Irregular Grain 

Let us consider the general case where the dipole is fixed in the grain body, and given by 

fi = ptiai + 112^.2 + /i3a3- (A28) 

where /ij are components of electric dipoles along three principal axes. In Paper I we disregarded the third term in 
Equation (A28) because of grain's axisymmetry. 
In the inertial coordinate system xyz (see Fig. 1), ai,a2 and as are described as 

ai = sin (f) sin (?x — cos cf) sin 6y + cos 0z, (A29) 
a2 = (cos (/) cos f/' — sin0sin^ cos0)x + (sin cos ^/^ + cos0sin^cos6')y + suiipsmOi, (A30) 
as = — (cos sin ^ + sin cos ^ cos ^)x + (— sin <;/) sin ^ + cos (j> cos ip cos ^)y + cos ^ sin 62,, (A31) 
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where (f>, and 6 arc Eulcr angles. 

Complex motion of the grain principal axes with respect to J results in an acceleration for dipole moment in the 
inertial coordinate system xyz: 



/i = /Uiai + /U2a2 + usais, 



(A32) 



where ai and a2 are given by 



ai 



^2 



+ 



+ 
+ 



as 



+ 
+ 
+ 



(^^ + r) sin (j)sm6 + 6^ cos (f) cos 6 

4>^ sin 6 cos + v cos (j)s\n.6 + 64> sin (j) cos 6 y — 9^ cos 6z, 
{(j)^ + tp'^ ) {cos (j) cos tp — sin sin cos 61) — 2(j)'tp{— sinc^sin-!/! + cos </) cos cos 
6^ sin (j) sin ip cos + 9(j) cos (j) sin sin 6 + 9tp sin ^ cos ^ sin 

— {<p + V'^ ) (sin cos tp + cos sin tp cos 9) — 2^tp (cos sin + sin cos cos 9) 
—9^ cos sin tp cos 9 + 9^ sin sin V' sin ^ — cos (p cos sin ^ 

— {tp"^ + u) smtp s\n9 + 9tp cos tp cos 9 

)^ + tp'^ ) {cos (p sin tp + sin cos -0 cos 6') + 2(ptp{s'm(pcostp + cos </) sin i/; cos 0) x 
9^ sin ^ cos ^ cos 9 + 9(p cos cos V' sin ^ — 9t^ sin ^ sin tp sin ^ 
(^^ + tj)^ ) (sin (psintp — cos ^ cos cos 9) — 2^tp (cos ^ cos — sin sin tp cos ^) 
— cos (p cos '0 cos 9 + 9<p sin cos V' sin ^ + cos (p sin V' sin 9 

— (-0^ + cos V' sin (9 — Oip sin i/i cos 6'j z, 



(A33) 



(A34) 



(A35) 



The precession and rotation rates ^ and tp are related to the angular velocity components as follows (Landau & 
Lifshitz 1976): 



By solving equations, we obtain 



ui = <f) cos 9 + tp, 
UJ2 = (pshi9 sintp + Ocostp, 
ujz = (p sin 9 cos tp — 9smtp. 

= LU2 costjj — LOg, sintp, 

U)2 sin tp + UJ3 cos tp _ jl20J2 + ^3^1 



sin^ 



7-2, ,2 , 7-2, ,2 ' 
tp = OJi — <p COS ^, 



(A36) 
(A37) 
(A38) 

(A39) 
(A40) 
(A41) 



Plugging into Equation (A32) with the usage of Equations (A33) and (A34) we obtain the acceleration components 
as functions of time. Performing Fourier transform for these components gives us the spectrum of electric dipole 
emission (see Fig. 2). 

The dipole emission power of this torque-free rotating grain can be obtained by averaging the /i^ over time: 



2 , 2^ 1 



3c3 



(A42) 



B. ELECTRIC DIPOLE DAMPING FOR DISK-LIKE GRAIN 



In the grain body system aia2a3, the dipole moment is given by Equation (A28). The orientation of axes aj in an 
inertial system are determined by Euler angles (see Fig. 1). The increase of grain's angular momentum over time 
arising from the acceleration of dipole emission is then 
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Using Equations (A28) and (A32) for (Bl) and averaging it over 4> and ip from to 27r due to torque- free motion, 
the non- vanished component is given by 

2 J3 

'3^/3 



dJz 
dt 



{cos"* - 3/i + 2) + cos^ Q{2,h - 2}^) + hf"} + ix{}f sin^ Q 



(B2) 



where we assumed /i^ = n\ = The components dJ.j./dt and dJy/dt arc averaged out to zero. 

For case 1 with /j.^^ : /i^ = 2 : 3, i.e., /x^ = /Xj = = /x^/3, we obtain (similar to HDLIO): 



dt 



2/^2 / 



[cos^ e{h^ - 3/i + 2) + cos^ e{-2h^ + 3h) + h^{l + sin^ 6)] 



For case 2 with /xy : /xj, = : 1, i.e., /xi = and /x^ = = /x2/2, 

= [cos" 9{h^ - 3/1 + 2) + cos^ e{-2h^ + 3h) + 

In dimensionless variables, we have 



dJ^ 
dt' 



2 J. 



'3 



3t' ' 



where r^^ = Ted,eff/TH,|| with 



Ted.eff = Ted,|| 



cos4 e{h^ - 3/i + 2) + cos2 6»(-2/i3 + 3/i) + + sin^ 6) ' 



3/,fc' 



2^3 



Ted.ll 



4A;BTgasM2 ' 



for case 1, and 



Ted,efI = Ted,|| X 



cos2 e{-2h^ + ih) + cos* 6»(/i3 _ 3/1 + 2) + /i^ ' 



Ted, II 



2fcBTgasM' ' 



(B3) 

(B4) 

(B5) 

(B6) 
(B7) 

(B8) 
(B9) 



for case 2. 

For ^ = 0, then Ted.eff = Ted,| 
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